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We study in this paper the ground state energy of a free bosonic theory on a finite interval of length 
R with either a pair of sine-Gordon type or a pair of Kondo type interactions at each boundary. This 
problem has potential applications in condensed matter (current through superconductor-Luttinger 
liquid-superconductor junctions) as well as in open string theory (tachyon condensation). While 
the application of Bethe ansatz techniques to this problem is in principle well known, considerable 
technical difficulties are encountered. These difficulties arise mainly from the way the bare couplings 
are encoded in the reflection matrices, and require complex analytic continuations, which we carry 
out in detail in a few cases. 



I. INTRODUCTION 



This paper is concerned with the study of 1 + 1 quantum field theories defined on a finite interval of length R. 
When the bulk is massless and conformal boundary conditions are chosen, such theories can be entirely understood 
using the formalism of boundary conformal field theory [1]. The case where one deviates from this situation either by 
adding an interaction in the bulk or at the boundaries (or both) is more difficult, and quite interesting as it involves 
crossovers depending on the bulk and boundary mass scales as well as the finite size R. 

Since the "conformal revolution", this topic was probably first considered in print in [2] (although it has been 
studied in unpublished work of A. Zamolodchikov) where the basic thermodynamic Bethe ansatz and Destri de Vega 
(DDV) approaches were delineated in the case of simple integrable theories, mostly massive in the bulk, with simple 
boundary conditions. Since then, the topic has attracted a fair amount of attention (see for instance [3], [4], [5], [6], 

In the present paper, we discuss in details the limit complementary to [2] where the bulk is a massless theory - 
mostly a free boson - and the perturbations sit on the boundary. We shall have mostly two theories in mind. The 
first is of the boundary sine-Gordon model type 1 
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In this last equation, refer to U q sl(2) spins in representations of spin ji,j r and q = e m P , The case with 

ji = j r = 1/2 could be called a "double Kondo problem". In that case, as usual, the /3 dependence of the interaction 
can be traded for a coupling anisotropy in spin space - notice that we restrict to the case where the I and r sides see 
the same anisotropy, a condition necessary to preserve integrability. 

This double Kondo theory musn't be confused with the two impurtity Kondo model. Although there also one would 
expect results to depend on the couplings for each impurity and the distance between the impurities, the detailed 
set-up is entirely different [8]. 

There are plenty of physical motivations to tackle these problems. On the condensed matter side, they can be 
thought of as generalizations of quantum impurity problems, where the bulk of the theory gets modified by the 
addition of a single local defect, e.g. a boundary or a quantum degree of freedom. Now, instead of only one impurity, 
we have two different ones coupled by a finite-size quantum field theory. One could label such a situation as a 
"quantum bridge problem" . An experimental setup for such problems could involve an interacting device - maybe 



1 In [2] the results referred to as boundary sine-Gordon model concern mostly a bulk sine-Gordon theory with Dirichlet type 
boundary conditions in 0, R. The terminology has somewhat changed since then, and boundary sine-Gordon theory nowadays 
refers rather to a free boson with a cosine perturbation at the boundary. 
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a carbon nanotubc or quantum wire - connected to two leads on either end. If the excitations of the device suffer 
from a gap when trying to infiltrate the leads, then in a first approximation all the action takes place at the contacts 
between the device and the leads. The excitations, upon hitting the contacts, will be reflected and mixed by the 
"boundary" according to some rules dictated by the specifics of the couplings, in other words the excitations will face 
some nontrivial, energy dependent boundary conditions, and the situation could be described by hamiltonians such 
as (1,2) in the low energy, universal limit. 

Possibly the simplest example of such a situation occurs when one considers two different (s-wave) superconductors 
connected by a bridge made of a quantum wire, e.g. a carbon nanotube, a construction called a Josephson junction. 
Since both superconductors inherently possess a gap for electronic excitations within their bulk, an effective theory 
can be obtained by considering only the wire itself with two "impurities" on cither side replacing the superconductors. 
Solving this effective model - which is essentially described by (1) leads to the determination of the Josephson current 
[9]- 

We do not know of any experimental condensed matter situation where the two boundary Kondo model would be 
relevant (although this does not seem impossible to imagine) . The model has however appeared in string theory in the 
study of tachyon instabilities [10] where the ends of the open string are coupled to a background, and the boundary 
spin acts in the space of Chan-Paton factors (and the free boson corresponds to the X 25 coordinate). The double 
boundary sine-Gordon model could presumably also be interpreted in the context of tachyon condensation following 
[11]. 

Notice that in the hamiltonian (2) we have not introduced the phase difference \- This is because such a phase can 
always be absorbed into a (gauge) redefinition of the spin operators S^. Alternatively, observe that a perturbative 
computation of the ground state energy of this model will only involve configurations of charges (the charges in the 
boundary vertex operators) which are independently neutral on the left and right sides, which leads to the possibility 
of shifting by a constant independently on each boundary (notice also that this ground state energy will then expand 
in powers of Af, A 2 .. The situation for (1) is very different, as all charge configurations which are only neutral overall 
(ie by combining left and right sectors) do contribute. The dependence of the ground state energy on \ 1S actually 
one of the main concerns of this paper. 

Let us now give some generalities about our approach. We consider a strip of width R and length L — > oo, with 
Euclidean coordinates x <E [0, R] and y e M. We will mostly discuss the double boundary sine-Gordon model in this 
paper (1), which belongs to a more general class of bulk and boundary perturbed conformal field theories, whose 
Euclidian action can be schematically written 

S = Scft+cbc + J dx J dy <S>(x,y)+ J dy($ Bl (y) + <P Br (y)) . (3) 

Here, we have as a starting point a CFT with action Soft, perturbed by some relevant operator <3>(x,y) contained 
within the operator content of that particular CFT. In the presence of boundaries, one needs to specify conformal 
boundary conditions (CBCs) at x = 0, R, and include the nontrivial boundary effects as a perturbation on the CBCs 
by some relevant boundary operators (in which we have included coupling constants gi, g r ) &Bi, r (y) = gi,r${%,y)\x=o,R 
at the left and right boundaries. 

The theory (1) turns out to be integrable. This can be shown by a simple generalization of the argument of [12], 
which we discuss here for completeness. The Lagrangian or perturbed CFT approaches both provide a symmetric 
stress-energy tensor T^ v with components T zz = T, T zz = satisfying the continuity equations 

d- z T = z e, z f = d- z Q. (4) 

In the presence of boundaries, as was shown by Cardy [1] , the choice of conformal boundary conditions ensures that 
no momentum can flow through the boundaries, in other words that the off-diagonal components of the stress-energy 
tensor vanish at the boundaries, T xy \ x= o,R = 0. In the presence of boundary perturbations, however, this equation 
gets modified to 

T xy \ x=0 ,R = -i(T - T)\ x =o,r = T^fyr (5) 

where i? is some local boundary field, which can be related to the boundary fields &b, r - 
Consider now the integrals over a closed contour C 

MQ= I (Tdz + Odz), P{C)= ( (Tdz + Odz). (6) 



c 



By the continuity equations (4), these do not change under deformations of the contour C: Pi(C) = Pi(C) = 0. Let 
us now take C to be the contour illustrated in figure 1. 
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FIG. 1. Contour of integration for the construction of the integrals of motion. 

We then can write 

o = Pi(C) + Pi(cr)= ]T p^c^ + P^Ca). 

a—u,d,l,r 

The contours C; and C r , by equation (5), can be seen to give the contributions 

Pl(Cl, r ) + Pl(C,, r ) = <V(t/l) - lVO/2). 

If we take y to be the time direction of our quantization scheme, this then shows that the quantity 
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is y-independent, in other words is an integral of motion. 

If the theory is intcgrable, these equations are but the first elements of an infinite series. Then, there are an infinite 
number of local fields T s+ i, s _i obeying 



d- z T s+1 =d x Q,- U 
in the bulk. Provided they obey the boundary conditions 



d- z f s+1 = c9 z e s _! 
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for some infinite subset {sb} € {s}, the theory will remain integrable despite the boundaries, with integrals of motion 

H ir(y) = / dx [ T s+i + ®s-i + f s+i + e.-i] + 0. ; i(y) + #s-Av)- ( 12 ) 

Jo 

Equations (11) define the possible integrable boundary conditions. They formed the fundamental set of conditions 
used in the analysis of [12]. 

Integrability can be used in two different - and equivalent through crossing - ways to calculate the ground state 
energy. To see this, let us switch to a hamiltonian formalism; imaginary time can run along x or y. For time flowing 
along y, it is natural to identify the first element of the infinite sequence of conserved quantities, equation (9), with 
the Hamiltonian. The boundaries are thus true boundaries in space, and the bulk excitations satisfy appropriate 
boundary conditions. In this language, baptized the L-channel, the partition function is obtained by tracing over all 
such states, with L — > 00 as the inverse temperature: 



Z = Tr e 



-LH, 



(13) 



On the other hand, in the i?-channel, we take time to flow along x, from an initial state to a final one (i.e. from x = R 
to x = 0). The Hamiltonian in the bulk is thus the usual one (without boundary contributions), and the associated 
space of states is the same as in the absence of boundaries. The partition function now becomes a simple matrix 
element, obtained by sandwiching the time-evolution operator H along time R between the initial and final states 
(the "boundary" states) \Bi, r ), i.e. 
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Z = (B l \e- HM \B r ), (14) 

This latter point of view leads to the TBA in the crossed channel, which we will use extensively in what follows. 
To handle (14) one uses the description of the integrable theory as a massive (if there is also a bulk perturbation) or 
massless (if the perturbation affects only the boundary) scattering theory. Such a theory is described by the spectrum 
of particles together with the bulk and boundary scattering matrices, and the constraint of factorizcd scattering. 

Consider for instance the massive case, and parametrize energy and momentum of particles through a rapidty 6. 
The explicit form of the boundary states is [12] 

\B) = exp [j™ deK a \6)A\{-e)A\{6)} |0) (15) 

where the amplitudes appearing in the exponential are given by the analytic continuation of the reflection matrices: 

K a \6) = R b J%-e). (16) 



Using this basic ingredient, a straightforward method to compute the ground state energy was proposed in [2], in 
the case where bulk and boundary scattering are diagonal. The non diagonal cases require more work, and discussing 
them is part of our purpose here. It so turned out however that even in the diagonal case, considerable surprises are 
encountered - as was first discovered in [4] - due in part to the way the bare couplings are encoded in the boundary 
scattering description. Calculating the correct ground state energy requires analytical continuation of the results of 
[2], a process which is quite difficult. 

The case of double Kondo model (2) is a bit different, as there are additional boundary degrees of freedom. The 
proof of integrability would proceed along similar lines, with appropriate modifications as discussed in [13]. The 
boundary state formalism can also be generalized to this case. We will, however, take a different approach based on 
the Destri-de Vega formalism. 

Most of this paper is devoted to the double boundary sine-Gordon model (sections 2 to 6). We will also present 
results for the double Kondo model in section 7. We shall first discuss analytic continuation by considering the case 
of the Ising model in section 2. In sections 3 and 4, we extend the analysis of [2] to the case where the boundary 
scattering is non diagonal while the bulk scattering is still diagonal. The analytical continuation is carried out as 
a simple generalization of section 2 for the case (3 2 — An (many of the results of this section have been previously 
published in [9], but we include a more thorough explanation here, which helps to understand the generalizations we 
propose). The other "reflectionless cases" require considerably more work; only the case f3 2 — is fully discussed, 
in section 5. The double Kondo model is tackled in section 6; there, the bulk and boundary interactions arc non 
diagonal, and we use a different method, based on the often called Destri-de Vega approach [14] (although the method 
has been used before in slightly different settings [15]). 

II. A SCALAR THEORY: THE ISING MODEL 

Let us now turn to a specific example, and consider the Ising model. As is well-known, the bulk action of this 
theory takes the form of a free action for Majorana fermions ip,^. At the left and right boundaries x — 0, R, there 
are two (in general different) values of the boundary magnetic field hi_ r coupled to the local boundary spin operators 

$l,r = ^{^P + "4>)x=0,R ai,r(v)- ( 17 ) 

where a; jr are boundary fermions, and we have chosen a convenient normalization. The total action is 

rR 



S = J dx J dy " i>d z $] + J dy(Ci + C r ), 



£l,r = ^i>i>\x=o,R + ai^ r d^ r + ht, r $i tr . (18) 

To proceed, let us do some elementary perturbative calculations, and adopt for a moment a more general setting. 
Consider thus a theory that is conformal in the bulk, and is subject to boundary perturbations at the left and right 
edges with coupling constants gi ir . We can write the partition function of our system as 



Z = (0\T V exp f dy [ 9l <f>(0, y) + g r $(R, y)] |0) 
Jo 



(19) 
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where |0) is the ground-state of the unperturbed system, T y is the y-ordering operator, and all the correlation functions 
are evaluated in the unperturbed CFT with the chosen CBCs. We introduce a conformal mapping by defining the 
new coordinates 

w = Wl +iw 2 = e i7Tz/R (20) 
in such a way that the original strip maps to the half-disc, as illustrated in figure 2. 




w i r PI 



FIG. 2. The conformal map between the original strip of u and the complex plane of z. 

The boundary is the horizontal axis, the perturbation for x = now sitting on the boundary w\ > 0, and the one 
for x = R sitting on the one for w\ < 0. Introducing d the dimension of the operators $ (such that the two point 
function on the boundary goes as l/|to| 2d ) and p = e~^ L / R , we find 

Z =<0|Tcxp^ 1 -^ (J^j [ 9l <S>(w)+g r 3>(-w)}\0> . (21) 



The strategy is to evaluate ItiZq for large L, and extract the leading term proportional to L, since in that limit one 
g inZo ~ _^; _ Elementary calculations give then 

2(1-0 
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This formula will be useful in what follows. 

Let us get back now to the particular case of the Ising model and consider the simplest case where one of the 
boundary fields vanishes, the other being equal to h. We then obtain 



where 



P 0lni /P N = 1 Jp<|xi|<...<|x 2 jv|<l X{ X 2 ' N 

in which < ... > c means that the correlators are evaluated in the unperturbed CFT. Evaluation to first two non trivial 
orders gives 

A8R 47r^n+l/2 + 8tt 2 ^ (n + l/2)(p+ l/2)(n + p+ 1) ( ' 

The first term is manifestly divergent (as would be seen as well from (22)), and requires some extra rcgularization. 
Since it diverges logarithmically, it is natural to expect any dimensionally regularized version of this term (such as 
given by the TBA) to go as ln(cst h 2 R). The next term, like all the subsequent ones, is well-defined. The sum can 
be evaluated in closed form: 

^ (n+l/2)(p+l/2)(n + p+l) = ?C(3) 
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giving rise to 
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Higher orders could be obtained without too much difficulty, but the foregoing expression will serve our purposes. 
Let us now compare it with the non perturbative results obtained from the TBA. The idea will be recalled in detail 
in the next section, so for now we just give the result, which can be found easily by taking the massless limit of [2]: 



k — h 2 R 
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It is not entirely straightforward to expand it in powers of h 2 R, in particular because the second order term actually 
involves a logarithmic piece. To proceed, if we write E = — j^I(x), x = h 2 R, we have 



AI = I{x) - 1(0) -- 
It is then more convenient to expand the derivative as 
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one finds then, after reintegration, that 

AI = x(c-l) + x\nx- J(-1)«M^ r dt( 
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These integrals cannot all be evaluated in closed form, but the lowest one can: 
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from which it follows that 
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This agrees with the perturbative calculation. In fact, one can also check the validity of the result to all orders (27) 
by comparing it with the lattice calculations in appendix B of [19], so things are quite satisfactory here. 

Surprises start when one considers the case of two non- vanishing applied fields hi , h r ^ . The TBA approach 
then gives 



E (hi,h r ) = - 



1 
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(30) 



This is confusing because it seems to have an expansion in powers of hf and h 2 only, a result in manifest contradiction 
with perturbative calculations. Related to this puzzle is the fact that (30) is the same for hih r > and hih r < 
while the two cases ought to be different based on even simpler arguments than perturbation theory: indeed, in the 
limit of large fields, the first case describes a partition function with fixed ++ boundary conditions, and E gs = —4^, 
while the one describes a partition function with fixed H — boundary conditions, and E gs = + the difference 

corresponding to the boundary dimension of the spin operator, d = ^. 

To make things more precise, the perturbative expansion (22) leads, at second order, to 



E (hi,h r ) = 
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h 2 h 2 1 
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(31) 
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To solve the apparent contradiction with (30) we attempt to expand it in powers of the magnetic fields. Since we 
already have a controlled expansion for one vanishing field, e.g. h r — 0, we consider 



E Q (hi,h r ) - Eo(hi,0) 



1 



This can be rewritten as 

E (hi,h r )- E (hi,0) 
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This is not analytical in h 2 ., as can easily be demonstrated by trying to calculate the derivative with respect to h~ at 
h r = 0, which diverges. To proceed, we make the change of variables k = h\h r Rx (assuming h\h r > 0),and get after 
some simple manipulation 
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Putting h r = in the bracket now gives a convergent integral, so 
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and thus, the first non trivial term in hih r is obtained by setting hi — in the integral 
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in agreement with the perturbative result. 

We thus resolve in a simple way the first paradox, by recognizing that the TBA integral must expand in fact powers 
of hi and h r , not hf, at least when both fields are positive (devising a scheme to perform this expansion is another 
matter) . 

There does remain a bigger problem, which is that the expansion of the TBA formula (30) is in powers of \hi\ and 
\h r \, and thus does not see the difference between the cases hih r > and hih r < (the cases hi, h r > and hi,h r < 
being meanwhile equivalent by spin reversal symmetry). 

The reason why we have this problem seems interesting. We do not think we understand it fully, but an important 
hint can be found by going for a while to the initial description based on the ordered, massive phase (that is, before 
we take the massless limit). The R matrix is calculated with a single (say, right) boundary at x = 0, but of course, 
the problem is undefined until one has set proper boundary conditions at (left) x = — oo. Since we are in the ordered 
phase, there are two possible such conditions, spin fixed up or down. For a given magnetic field on the (right) 
boundary, the two (left) boundary conditions at x = — oo should correspond to different problems, with presumably 
different states describing the right boundary. Ghoshal and Zamolodchikov in [12] argue that spin up and hi > 
correspond to the usual boundary state, while spin down and hi > correspond instead to an "excited boundary 
state", obtained formally by some procedure of analytic continuation. We believe that a similar feature holds, the 
other boundary in our finite geometry playing the role of the boundary condition at x = — oo in [12] . If this is the case, 
the TBA approach proposed in [2] lacks a crucial element, and is good only when, roughly speaking, the boundary 
couplings would not lead to existence of frustration in the integrable massive bulk theory compatible with the massless 
boundary perturbations. In the Ising case: when the two fields have the same sign. 

How to handle the case hih r < thus requires some additional ingredient. To see what it may be, let us calculate 
the ground state energy of the Ising model on the strip by brute force, diagonalizing the hamiltonian H. Following 
calculations initiated in [17] we have 



k<0 



The modes k are determined from the quantization equation 

2ik hi + 2ik 



hf 



2ik hl-2ik 



= -1. 



(36) 



(37) 



One can easily show that this equation has only real k solutions, and that if k is a solution, so is — k. Therefore, the 
sum in (36) runs only over the k > values. 
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Now a detailed study of the quantization equation shows that all solutions are analytical functions of hf and h% 
except for one of them, which in particular goes as k = ±fco w ±||/i;/i r | at small magnetic fields. 

Manipulations of (36) in the usual way lead to the TBA integral formula (30). Suppose now that we believe this 
result for hih r > 0. To get the result for hih r < 0, we could take the formal expansion of (30) in powers of hi, h r and 
continue analytically in the variable hih r . This, it seems from the previous analysis, involves only the modes ±fco- If 
the mode +ko is included in the sum for hih r > 0, the analytically continued result will be obtained by discarding it, 
and using instead the mode — fco: in other words, we suggest that the energy for negative hih r simply reads 

k — hfR k — tilR 

1 + 



1 f°° 

Eo(hi,hr) = J dn\n 



+ k , (38) 



K + hfR K + hfR 

If hih r < 0, this expression goes as —\\hih r \ + \ \hih r \ = \ \hih r \ = — jhih r , as desired. We also observe that at large 
magnetic fields, it goes to E « — 4^ + ^ as k — > in agreement with the conformal prediction for +— boundary 
conditions. 

Two tests are therefore satisfied. We also observe that (38) implies the simple result 

E (hi,hr)-E (hi,-hr) = -ko. (39) 

Note that the idea of adding an extra contribution in different regimes of couplings appears first in [4] . 

To provide supplementary justification for this construction, we have performed some numerical checks on equation 
(39) using Mathematica. The plots were produced using exact diagonalization of a lattice Ising model at its bulk 
critical point, for system sizes up to ten sites, with equal absolute values for the left and right boundary magnetic 
fields. Both cases hi — h r " to the root k . Figure 

3 shows our results. 




FIG. 3. Ground-state energy difference for the two-boundary Ising model with boundary magnetic fields hi — h r and 
hi = —h r . This is plotted as a function of the scaling parameter B}^ 2 h for system sizes up to ten lattice sites, and compared 
to the theoretical prediction in equation (39), i.e. the root ho- 



Difficulties related with the relative values of the two boundary parameters will be generic - and worse - in the 
case of the double boundary sine-Gordon model, and require considerable efforts. On top of this, an additional 
difficulty arises due to the presence of the phase \'- even if t ne bulk is chosen to be at the reflectionless points and 
hence have diagonal scattering, the boundary reflections cannot be simultaneously diagonalized, requiring substantial 
modification to the crossed channel TBA of [2]. We would like to illustrate this in the case of the free fcrmion point 
1 = 4tt. 



III. ALL THE DIFFICULTIES IN A NUTSHELL: THE DOUBLE BOUNDARY SINE-GORDON MODEL 

AT [3 2 = 4tt. 

A. The crossed channel TBA 

At the free fcrmion point 1 = 4ir, the (massive) sine-Gordon theory simplifies considerably. The only fundamental 
excitations are two fermions, the soliton and the antisoliton. Their scattering matrix in the bulk is trivial, independent 
of rapidity - the creation and annihilation operators obey the canonical anticommutation relations 
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{A a {e),A\{6')} = 6 ab 5{8-6') (40) 

with all other anticommutators vanishing, a = ± labelling soliton/antisoliton particle types. Scattering at the 
boundaries is encoded in the R matrix 



R++ R+-\ = P+ Q 
R-+ R— J - \ Q P- 



(41) 



There will be one reflection matrix for each boundary, and its elements will depend on the rapidity of the incident 
particle, as well as the bulk and boundary parameters of the field theory. It is easier to obtain the reflection amplitudes 
formally by solving the boundary bootstrap equations; one is thus left with some undetermined relation between the 
parameters entering the R matrix, and the parameters in the field theory. We will have to tackle this issue in more 
details later. For now, we observe (as discussed in details in the appendix) that the point 1 = 47r being equivalent 
to a free fermion theory, explicit calculations can easily be carried out, and one finds [22] 

Pi(9) = [coshS- in,rCoeh(«Ti*i,r)]/A,r(»), 

- = <" 2 > 

n ia\ ■ ,6-i4> Lr -iir/2 . 8 + i<f> Lr + in/2 2 

A,r(#) = -m,r cosh( - ) smh( - ) - cosh 9, (43) 

where 7; ir cx Af r /m and (f>i = 0, <f> r — 4>o- 

In the i?-channel description of the problem, the explicit form of the boundary states then reads 

\Bi, r >cx exp (jf° MKft(0)Ati-0)Al(0)) |0 > . (44) 

Here, the K matrices are given by the analytic continuation of the boundary R-matrix according to 

Ki AO) = 7^ s t r x 

sinh 2 6 + 7i , r cosh (a^rj cos h (t^hr^ 

sinh20 — 2 sinh# + sinh(6* — i<f>i ir ) \ 

-2smh6 + ji, r smh(8 + i(f)i, r ) sinh 29 )' ^ > 

As it should, this obeys the boundary cross-unitarity condition K ab {9) = —K ba (—9), in view of the fact that the 
scattering matrix is simply —1. 

It is now straightforward to take the massless limit through 

m^O, ^o^oo, me "=me W9 '=/, n finite. (46) 

8' (which we will often denote by 9 when no ambiguity is possible) is then integrated from — oo to oo. We will often 
denote fie 8 = k, and in practice set /i = 1 (sometimes the massless limit is taken with /j = 2 instead, which does not 
affect any of the results. Also, we use the notation \x briefly in appendix C, but no confusion should be possible). 
Performing the massless limit on the if- matrix (45), we get 

KiAO') = ,. ' A2 Lrl. i- 147, 




(Note that some care has to be exercised in getting the K matrix in the massless limit directly from the R matrix in 
the same limit). The couplings A are renormalized couplings, in this case A = i/ttK. 

Equipped with this data, we now consider the calculation of the ground state energy generalizing [2] . Although this 
was already carried out in [9], we include here some more extensive explanations, which will be useful in the following 
sections. 

Let us now turn to the calculation of the partition function in the limit L — > oo, where the boundary states are 
well-defined. Our starting point is a version of equation (14): 

z = y(B l \a)(a\B r ) e _ RE 
^ (a a) 

a 

where a is an index running through all states with nonzero inner product with the boundary states. In view of the 
structure of the boundary states, we will thus need to trace over all states of the form 
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\a)e{Aij-e N )4je N )...Al(-e 1 )Al i (e 1 )\o), va e n,v9 



(49) 



We have to be very careful here: in the case of a theory with only one type of particle, the rapidities can be ordered 
using the strict inequality On > ... > 6\, since the Bethe wave functions vanish identically for coinciding rapidities 
(see e.g. [2]). Now, however, we have more than one species of particles, so we can have coinciding rapidities without 
having vanishing Bethe wave functions, as long as the labels of particles at coinciding rapidities do not get reproduced. 
A more economical strategy consists in introducing creation operators for particle pairs, and a quartet one [9]: 



PLW = Ai(-6)Al(0), Qt(0) ee A\(-9)A\(9)A\{-9)A\{0), 9 > 



(50) 



where the pair creation operators come in four flavours, labeled 11,22,12 and 21. We can now consistently restrict 
the rapidities to be positive and noncoincident in the set of states \a). Note that because of basic anticommutation 
relations and the Pauli principle, the structure of the boundary state dictates that the quartet appears with fugacity 
K 11 {9)K 22 {9) - K 12 (9)K 21 (9) =detK. 

We can now rewrite our trace over states in (48) as running over the set of all states of the form 



l«) e { Pl Nb J0 N )...Pl bi (9 1 )QH9%J...QH9l)\0), 
V A, A g e N, 9 N > ... > U 9% > ... > 6\, 9i ? 9 q A. 



(51) 



For a state with A pairs with rapidities 6i,...,9n and N q quartets with rapidities 6\, 6 N we can easily show by 
explicitly expanding the exponentials in the expressions for the boundary states that 



{Bi\a){a\B r ) 
(a\a) 



N 



N 



= ]jK^(9 l )K^(9 l )]Jdet[K l (9])K r (9]) 



(52) 



To construct the entropy of a field configuration, necessary to complete the TBA, we use the following argument. 
Imagine putting our system in a very large but finite domain of length L in the space-like direction of the R-channcl 
picture. Since all bulk scattering is given by the trivial scattering matrix S = — 1, the quantization rules dictate that 
the allowed rapidities appearing in the pair and quartet creation operators lie on a lattice of almost evenly spaced 
values (separated by intervals of order 1/L, and symmetric with respect to sign exchange). Moreover, in view of the 
Pauli principle, no superpositions of creation operators can take place on a given rapidity. Thus, the entropy of a 
state configuration with N a b pairs of type a, b and of N q quartets is very simple to calculate, being that of a gas of 
mutually excluding particles of five different types. Explicitly, for At total allowed rapidities, we have the simple 
formula 



N 11 \N22lN 12 \N2ilN q \(N T - An - A 22 - A 12 - N 2 i - N q )\ 



Introducing the usual occupation densities p, p a b, p q as well as the hole density ph = p 
limitc L — > oo, we can rewrite this in the Stirling approximation as 



S^L I d9 



■f 

Jo 



p\np - ^Pab In Pab - Pq hip, - Ph lnph 



ab 



(53) 



Y^ab - Pq as we take the 



(54) 



The energy is on the other hand given by 



E 



Jo 



d9 



+ 2Pq 



ab 



2m cosh#, 



allowing us to put everything together to obtain the partition function 

Z w J [dp] expj^^ de (lnKf a Kf - 2mi?cosh6>) p ab + (]ndeb[K t K r ] - 4mRcosh9) p q + S[p] 



(55) 



(56) 



Introducing the quasienergies e _£ab = p ab /Ph and a similar one for the quartets, we can obtain the TBA equations 
by using a saddle-point approximation on the partition function. From the quantization equation we get the con- 
straint ph{\ + ^2 ab e~ eab + e~ £ «) = ^ cosh#, which when substituted in the saddle-point equation leads (after some 
straightforward manipulations) to 



10 



InZ 



^ [ d6 cosh0ln( l + Ve-^+ e - e O 



e a6 = 2m.Rcosh0 - In K^K™, 
e q = AmR cosh 6 - \ndct{KiK r ). 

Reinterpreting this as the ground-state energy of the original problem, we arrive at the final formula 



En — 

2tt 



d9m cosh# In 



1 + triK^e 



-2mR cosh 



+ det{KiK r )e 



-AmR cosh 6 



where the K matrices are given in equation (45). 
We can now take the massless limit using the results (47) to obtain the final expression 



— r 



dn In 



i + 2 K 2 + AfA 2 r cos2x c _ KR + (k - Aj)(n - A 2 ) ^_ 2kR 



{n + Af)(n + AD 



( K + A?)( K + At) 



(57) 



(58) 



(59) 



In the first appendix, we provide the details of an independent derivation of this result using refermionization, 
confirming that the TBA construction that we have employed is justified. 



B. Analytic structure of the ground-state energy 



The study of the analytical properties of the ground-state energy of the double-boundary Ising model in the first 
part of this paper has shown how careful one has to be in interpreting the data coming from the TBA equations. It 
is here again very easy to miss out on the correct analytic continuation procedure to be imposed. 

Let us thus carefully study our expression for the ground-state energy of the double-boundary sine-Gordon model 
at the free fermion point, equation (59), whose bulk term we reproduce here: 



1 f°° 

E = / dn In 

4vr Jo 



K 2 + A 2 A 2 cos2 X R («-Af)(«-Ag) 2kR 
( K + Af)( K + A2) + ( K + Af)( K + A2) 



(60) 



First of all, we expect the ground state energy to be an even function of the phase Xi as 1S obvious from the 
boundary sine-Gordon Hamiltonian. We can consequently concentrate on positive \ from now on. The first crucial 
observation is that equation (60) is valid only in the domain \ € [0, 7r/2], since it exhibits a singularity at \ = 7r /2: 
for x € [0, 7r/2[, the argument of the logarithm has zeroes in the complex plane of the variable k that hit the real axis 
at k — when x reaches ir/2. 

We can gain a bit more intuition about this formula by considering the limits A/,A r — > oo, where the integral 
reduces to 

1 f°° 

; = / dxlnri + 2cos2xe- :E + e- 2:l; l . (61) 

47T-K Jo 



This integral is tabulated, and we get 



As a function of the phase difference \, the basic integral representation (59) for the ground-state energy thus has a 
cusp at x = tt/2. 

A simple computation shows that this should not be the case. Namely, let us consider the boundary Hamiltonian 
(A12) in the limit A/. r — > oo. In that case, the field becomes pinned at x = to <fi(0,t) = 2y / 7r(n + 1/2), while 
at x = R it gets pinned to 2y / 7r(m + 1/2) + xl\f^i where n,m e Z. This immediately tells us that the zero mode 
operator ITo has to have eigenvalue n = 2y / 7r(p — 2^)iP e ^- F° r tnc phase \ in the fundamental domain \\\ < it, 
the lowest-energy state is thus the one with p = 0, so the energy simply becomes (using (A6) and adding specifically 
the contribution from the dynamical modes) 

exactly. This thus clearly demonstrates that the basic integral representation (59) works in the domain \\\ < n/2, 
but fails for n/2 < \x\ < 7r. 
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Let us now study the case of small A^ r . From the integral representation for the ground-state energy (59), we can 
write (after factorizing (1 + e~ K ) 2 in the logarithm and a few more basic manipulations) 



-En — 



7T 

24R 



4nR 



dn In 



1 - 



4t](k(1 + e K ) + 77 sin 2 x)e 



(1 



1 (k + T]) 2 



(64) 



where we have for simplicity considered A; = A r and defined the dimensionless parameter 77 = Af r R. One sees that 
it is not possible to do a simple Taylor expansion in 7/ in this formula. Rather, to evaluate it in this limit, we observe 
that the main contributions come from the region of integration where k ~ 77. Defining £ such that 77 << ^ << 1, we 
can approximate the integral term as 



4nR 



dn In 



rj(2n + r\ sin 2 x) 



1 



4nR 



dn In 



r k 2 + rj 2 cos 2 % 



(K + 7?) 2 



(65) 



This can in turn be approximated by rescaling by 77 and replacing the upper bound on the integral by infinity. The 
resulting integral is tabulated, and we finally obtain 



A? 

l.r 



cosxl + est. 



(66) 



The absolute value gives rise to the same kind of cusp as the one obtained in the limit A; !r — > 00. Here, we can 
compare with straightforward perturbation theory which gives us cosx instead of its absolute value. 

These considerations again indicate the presence of a problem with the integral representation (60) around the point 
X = 7r/2. In order to repair this, we can first use a procedure very similar to the one used in the Ising case. Namely, 
we identify the ground-state energy Eq with the integral representation (60) only in the domain |x| < tt/2. The 
domain tt/2 < |x| < tt is accessed by analytic continuation, which is performed based on the explicit diagonalization 
(see appendix A). First, note that all the roots of the quantization condition (A25) are well-behaved as a function 
of x except for the pair closest to the origin, which we label as ±fcn- This is the root that we need to analytically 
continue, as in the Ising case. Thus, in the domain 7r/2 < |x| < 7r, the true ground-state energy should be given by 
the integral representation (60) plus fen. 

The validity of this procedure can again be checked in the limits of strong and weak boundary pairings. In the 
limit A;. r — ► 00, the root fcn is given by fen = x ^ 2 ■ This leads, in the domain n/2 < \\\ < n, to a ground-state 

energy of E = + — h x ~r 2 = 517? + 2ttK> wn i cn i s tnc result required by the above considerations using 
the zero-mode operators. In the weak boundary pairing limit Aj >r — > 0, the root fcn behaves as — ^A 2 r cosx, leading 



in the domain tt/2 < |x| < tt to Eq - 
perturbation theory. 

To summarize, we write our final result as 



A 2 

-r I cos xl 



|A 2 r cosx + est. 



cosx + cs t -j as required by simple 



£ SG (A ; ,A r ,x) =£ (A ; ,A r ,x), 
E^ G (A h A r , x) = E (A h A r ,x) + fcn, 



Ixl < tt/2, 

TT/2 < Ixl < 7T, 



(67) 



with 2tt periodicity in x> aud where Eq is the integral (60). Plots of this ground-state energy as a function of the 
phase difference x can be found in [9] for various values of the boundary parameters. 

To extend the method away from the free fermions point, it is necessary first to reinterpret the addition of fco 
directly within the TBA approach; this will help make the case A = 2 more transparent. We start by writing down 
again our results from the TBA. For A = 1, since the scattering in the bulk is trivial, the quasienergies obey the trivial 
equation 



(68) 



where the ground-state energy has the integral representation (we have rescaled the boundary parameters according 
to A LiR = R 1 / 2 Ai tr ) 



Eq 



— r 



dn In 



1 



; 2 + A|A|cos2x 



'(k + A 2 l )( k - 



6(«) 



(«■ 



Ai)(«- 



A|) 



(«■ 



A|)(k + A 2 ,) 



-2e(«) 



(69) 



and we have performed some simple rescaling. Now of course, the equation (68) defines the quasienergy everywhere 
in the complex plane of k. In particular, e(n) is pure imaginary if and only if k is also pure imaginary. Moreover, 
this equation does not need to be modified as a function of x, i-e. there is no analytical continuation needed for the 
quasienergy. 
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If, however, we closely examine the ground-state energy, it is clear that something goes wrong. Looking at the 
integrand at k — yields 



In 



l + 2cos2xe~ ew +e 



e(0) , „-2£(0) 



2 In cos x (70) 



which clearly diverges when \ — * 7r /2, so wc know that there is a singularity hitting n = at \ = ir/2. For general 
X, the position of this singularity is found by selecting the proper solution of the singularity condition 

2e((e) ^+A|A|cos2x (K) ( K -Aj)(«-A^) _ 

+ '( K + A|)( K + A| J ) e + ( k + A|)(k + A|) U - [n) 

A careful look at this shows that we must look for such a singularity at values of k for which e(«) is purely imaginary, 
and therefore, according to what we wrote before, on the imaginary axis of k. We therefore note the singularity 
position as in s (note that if ik s is a singularity, so is —in s ). The picture is then as follows: singularities at ±in s 
touch zero at \ = 7r /2, and the ground-state energy integral (69) must be modified to take this into account. Simple 
contour-integral manipulations then allow to write the modification to the ground-state energy, which should be added 
for x > 7r/2, in the form SE — ^ making contact with the previous mode expansion arguments. 

This is how the analytical continuation should be understood for the case A = 1 in the language of the TBA, and the 
logic of this procedure will allow us to find the solution to the analytic continuation problem in the case A = 2. There, 
things will be more complicated: in particular, the self-consistency equations will themselves have to be analytically 
continued, and the singularity positions will lie in a different domain of the complex plane. 

IV. THE DOUBLE BOUNDARY SINE-GORDON MODEL: SOME GENERAL RESULTS AT THE 

REFLECTIONLESS POINTS. 

A. The ground state energy 

The massive sine-Gordon theory has its discrete symmetry 0^0+^/3, iVeN spontaneously broken for 1 < 8n. 
The spectrum has a number of fundamental particles classified as a soliton-antisoliton pair, and breathers (bound 
states). Although the soliton and the antisoliton are always present, the number of breathers is dictated by the value 
of the interaction parameter (5. Namely, the breathers are labeled with the integer n s.t. n = 1, 2, ... < A where 

A=^"l- (72) 

Solitons and antisolitons carry respectively positive and negative topological charge, whereas the breathers are neutral. 
If the soliton mass is given by m, then the breathers carry masses 

T17T 

m n = 2msin — — , n = 1, 2, ... < A. (73) 
2\ 

Let us adopt like in the free fermion case the notation (0) and Al_ (0) for the creation operators at rapidity 6 of 
the soliton and antisoliton, respectively. The creation operators for the n-th breather at rapidity will be denoted 
B\{6). Factorized scattering in terms of these is described by the fundamental scattering amplitudes 

4(^)4(02) = O (0i - 2 )Al(0 2 )Al(0 1 ), 

A\_(6i)A[_ (62) = b{0 1 - 2 )Al{0 2 )Al(0 l ) + 0(0, - 9 2 )A\(6 2 )Al(e{). (74) 

The amplitudes for these scattering processes are well-known [21], and take the form either of infinite products of 
gamma functions or of more compact integral representations. These are reproduced in appendix B. 

The introduction of a boundary interaction is described like in the free fermion case by a reflection matrix on the 
left and right boundaries respectively. In addition, there is an amplitude B n (0) for each breather - the scattering 
cannot mix particles, except solitons and antisolitons. 

The boundary boostrap method applied to this case yields the four amplitudes P±, Q± in terms of two free param- 
eters £, k which we will call "IR" parameters following [6]. Their relationship to the "UV" parameters A, \ appearing 
in the boundary Lagrangian is a complicated issue in general, which we will address later on. Suffice to say now 
that the results of the boundary bootstrap for the boundary scattering amplitudes are summarized in appendix C. In 
particular, we have (3(/)o/2 = X- 

We will restrict in what follows to rcflcctionless points where the bulk scattering is diagonal (that is, c(0) = 0). They 
occur at A e N. In these cases, there will be A — 1 breathers in the theory, with the original interaction parameter given 
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by P 2 = • As in the free fermion case, the solitons and antisolitons have to be carefully treated when ordering the 
rapidities. As detailed before, it is better to move to a basis where only positive rapidites are used, and the states 
involved in the trace for the partition function are constructed using pair and quartet creation operators P\ b , (see 
the discussion around equation (50)). Here, we will need to add the breather pair creation operator, defined as 



4(0) = i?t(-0)I?t(0), e>0. 
Now, we simply need to add all the possible breather pair states to our set of states to trace over: 

l«> e { ^ Nfew (^)..< fel ^i)Q t (^)...Q t W)^ Ni) (^)"-<(^)l )' 
v N,N q ,N b e n, e N > ... > e u e q Nq > ... > e\, e t ^ o]}. 



(75) 



(76) 



In order to compute the partition function in the thermodynamic limit, let us start with the quantization conditions 
for the allowed rapidities of the different types of excitations. For definiteness, let us consider a state with N pairs 
with rapidities 9\,...,6n, N q quartets with rapidities 9\,...,9 q N and Nb breather pairs with rapidities 6\,...,6 Nb . 
Let us imagine moving the particle in pair i having rapidity 9i around the system, and requiring periodicity of the 
wavefunction. This leads to the quantization condition 



N 



N„ 



e imLsinhei S(26i) 11 s{e i -e j )S{6 i + e j )J{s 2 {e i -ei)s\6 i + ei)J{s^\e i -e b l )s^\e i + e b l ) = i. (77) 



fe=i 



z=i 



The other quantization condition is obtained by moving one member of the breather pairs around the system in a 
similar fashion: 



N 



Am n sinh 9' 



5-K.«i)(26i 8 b ) Y[ S^(6 b i - 9j)S( ni \0 b + Oj) x 



N„ 



Y[ s* ( ™ i) {e\ - e q k )s (7M) (e\ + &i) Jl s^ ni ' ni \e b i - e b )s {n ^ ni) {e\ + e\) = 1. 

i&i)=i 



(78) 



k=l 



As usual, we define densitites of occupied and unoccupied rapidities but now for pairs, quartets and breather pairs. 
A complete set of states is obtained by considering positive rapidities only (remember: our pairs have one component 
at positive rapidity, and one at negative rapidity). For example, [p a i, + p^ b ]Ld9 gives the number of allowed rapidities 
in the interval 9,9 + 69 for pairs of type ab, with p ab the density of filled states. The symmetry properties p a b(~9) = 
Pab(9), p q {—9) = p q {9) and p n (—9) = p n {6) can then be used to simplify the resulting equations. In logarithmic form, 
we get (using the usual notation a * b{9) = J_ d9'a(9 — 9')b{9')) 



Y,P"b(Q)+P q (0)+P h (°) = £cOsh#-$* 



ab 



P„(0) + ^(0) = ^cosh0 -*<»>* 



J2Pab(0) + 2p q (9) 



A-l 



71 = 1 



Y,Pab(9)+2p q (9) 



ab 



A-l 



where the integration kernels are obtained from the scattering matrices according to 



*'*"<"> = 27,15'" sM(l »- 



v ' 2nid9 



(0), 



(79) 



(80) 



with the explicit form of the scattering matrices to be found in appendix B. 

The boundary matrices will appear in an interesting fashion in the TBA as rapidity dependent fugacities. For 
example, for a state with N pairs with rapidities 9\,...,9n, N q quartets with rapidities 9f, ■■■,9 q N ^ and Nb breather 

pairs with rapidities 9\,...,9 b N , the fugacity is 



(£?,|q)(tt|B r ; 
(a\a) 



N Nq N b 

]jK b ^(9 l )K^(9 l )]Jdet[K l (9 q )^^^ 

3 = 1 



(81) 



1=1 



fc=l 
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In the thermodynamic limit, this expression can be rewritten in terms of densities as 



llL ^\a)(a\B r)=L 
(a\a) 



dO 



_ ab 



A-l 



+ In det [Kt (6)K r (%(«)+^ \n[K\ n) {6)K™ {6)] Pn (9) 
In a similar way, the energy can be written as 



n=l 



(82) 



E = L d9 



o 



A-l 



2mcosh6>(^ Pab (9) + 2 Pq (9)) +^2 2m n cosh 8 p n {6) 



ab n=l 

whereas the entropy becomes in the Stirling approximation as L — > oo 



(83) 



/ 

Jo 



S^L d9 



P ab + Pq + Pq + P k ) m (X! Pab + Pq + P k ) Pab 111 P ab 



ab 



ab 



ab 



A-l 



Pq In Pq - P k In p h + ^ [(Pn + Pn) HPn + Pn) - Pn In Pn ~ P„ In Pn \] 



n=l 



(84) 



Equations (82, 83, 84) can now be substituted in the expression for the partition function in the R-channel, (48), which 
now becomes a functional integral over the various densities. The TBA equations then arise from the saddle-point 
evaluation of this expression, subject to the constraints coming from the quantization conditions (79). Introducing 
the usual quasienergies, we obtain after standard manipulations 



€ ab {6) = 2mRcosh6 - \nK b l a {9)Kf{9) + $ * In 
e q (6) = AmRcoshO - lndctKi(9)K r {9) + 2$ * In 



cd 



A-l 



1 + + e -««W 



cd 



e n (9) = 2m n Rcosh0 - \nK[ n) {6)K^{6) + $ (n) * In 



+ ^*( n )*ln[l + e - £ »W], 

n=l 

A-l 

+ 2^$(") *ln[l + e - £ »W], 

71=1 

A-l 

^ $(»."») *ln[l + e - £ »W]. (85) 



m=l 



Substituting these back in the expression for the partition function, and absorbing the boundary factors into a 
redefinition of e's, we arrive at the ground-state energy 



1 f°° 

E = / d9mcosh9ln 

2tt J 



1 + trKi(e)K r (9)e- e ^ + det Ki{9)K r (e)e- 2eW 



1 f 00 r 

/ d0Vm„cosh01n 1 + K\ n) {fi)K^ {9)e~^ 

2nJ ° ti 1 



(86) 



where the energies are self-consistent solutions to the system of equations 

e(9) = 2mRcosh9 + $ * In [l + trKi{9)K r {9)e- e{e) + det Ki(9)K r (9)e' 2e{e) 



+ 



A-l 



■ ^ $ (n) * In 1 + K\ n) (9)K^ (9)e 



71=1 



e n {9) = 2m„i?cosh<9 + $ (n) * In [l + trKi(6)K r (9)e- eW + det Ki{6)K r {9)t 



-2e(9) 



+ 



A-l 



+ $ ( "' m) * In 1 + Ki m \9)K^(9)e 



m—l 



71=1, A — 1. 



(87) 
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B. Universal form 



These equations can be put in a simpler, universal form, following [29]. The advantage of this reformulation is that 
only one integral kernel remains. For a given value of A, we label the breathers by 6=1,..., A — 1, the soliton by 
b = A, and the anti-soliton by b = A + 1. We can then use the remarkable identity [29] 



[Sab + ®ab(k)} 1 = 5 ab 



1 



2 coshfc//i 



Lab 



(88) 



holding between the Fourier transform of the scattering kernels and the incidence matrix L ao of the relevant Dynkin 
diagram (for sine-Gordon at A, this is the diagram of the D\ + \ algebra). The Coxeter number is h = 2A in that case. 
After basic manipulations, we obtain the universal form 



e a (0) = v a (0) + J2 L ^u * [-VbW + *b(0) + lnF 6 (0)] , 



(89) 



where the universal kernel is 



A 



2tt cosh \9 

and the boundary dependence appears in the F functions 



Fb(9) = 



l + K\ h) {6)K?\e)e-^ e \ 



[1 + trKi(0)K r (0)e- e W + detKi(8)K r (9)e- 2 < e 1] 1/2 , b = A, A + 1. 



b = 1,...,A-1, 



(90) 



(91) 



in which we have defined the kink and antikink quasienergies as e\ — €a+i = e. Finally, the driving terms v are 
redundant in the equations (the v a gets cancelled by the ^2 b L a b$ u * [— vt\), but very useful in the numerics, where 
they provide the correct asymptotics at large rapiditcs. They are 



v a (9) = 2m a Rcosh9. 



(92) 



C. The massless limit 



Taking the massless limit as described earlier leads to the following equations: 
1 f°° r 

E = / dKlnll + trKii^Kri^e-^ +detKi(K)K r ( K )e- 2£f(K 

4tt Jo 
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where the functions e, e„ are self-consistent solutions to the system of equations 
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n = 1, A — 1. 



(93) 



(94) 



The integral kernels are obtained from the massless limit of the known ones, and are given explicitly by 
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The boundary scattering matrices become in the massless limit 
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(95) 
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We have found it convenient to use a renormalizaed parameter A in this expression. What we call A A+1 coincides 
with what is called in other works [18] Tg. This can be checked by comparing with the R matrix in the massless 
sine-Gordon model studied in [18], where the only difference with the present case is that the phase is missing. Now 
one knows from elsewhere the correspondence between the bare parameters A and Tb'- 



„ ^ A±i 
T R = A— 



T(1/2A) 



We thus get the correspondence 

A = A 



r(i/2)r(i/2 + i/2A) 
r(i/2A) 



sin 



A + l 
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1/A+l 



sin 



A + l 



r(A/(A + i)) 



r(A/(A + i)) 
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(97) 



(98) 



This reduces to A = y/irA in the case A = 1, (i 2 = A-k as already discussed. 

The additional piece of information that is not in [18] is the relation between the phase in the K matrix and the 
phase in the hamiltonian: one finds (see the apendix) r\ = (A + 

The trace and det of left and right K matrices then read explicitly (we have used iji = 0, r] r = (A + l)x) 



trKi{K)K r {n) = 2 



detKi(K)K r (n) 



The breather K matrices are 



e 2A 
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^){n + A r x e l ~> 



[k2 a + a a+i a a+i cos(a + 1)x]; 

-2 
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~n 2X + (-1) A A 2(A+1) 1 \n 2X + (-l) x A 2 J x +V 



(99) 
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(100) 



Notice that the phase \ appears only in the K matrices for solitons and antisolitons, and for a single boundary, could 
easily be gauged away by a redefinition of the soliton and antisoliton. 
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V. THE DOUBLE BOUNDARY SINE-GORDON MODEL FOR A = 2 



A. Generalities 



Let us now turn to a precise investigation of what the TBA equations entail. For the particular case A = 2 
(corresponding to j3 2 = 87r/3, which is the next simplest case to tackle after the free fermion point A = 1, 1 = An), 
we can use the preceding results to rewrite the TBA in a more compact form. We have explicitly (after rescaling to 
the more convenient dimensionless parameters RA^ 2 = A"^) 



i z^oo fK roo 

/ dn\n[l + W 2 {K)] - -1— / d/sln[l + B 2 (k) 
Jo A7tK Jo 



W 2 (k) 



4nR 

2(k 4 + A|A3jCos3x)e- e2 ( K ) + (k 2 + A| - V2kA 3 l /2 )(k 2 + A| - ^kA^ 2 ^ 2 ^ 



(k 2 + A3 + V2kA^ /2 )(k 2 + A 3 R + V2nA R /2 ) 



_ A 3/2 _ A 3/2 

Bo(k) = - f=fl_ e -«w fion 

2i j K + Af« + Af ( J 

with the quasienergies determined from the universal form of the massless TBA equations 

o r°° K * K 2 

ei(«) = - / d/c'-J {2e 2 («') + ln[l + W 2 (k')}]} = U * {2e 2 ( K ) + ln[l + W 2 (k)}} , 

7T Jo K + K 

e 2 («) = - / d/c'-j { ei («') + ln[l + S 2 (/t')]] ee U ? [ Ci (k) + In[l + B 2 (k)]} , (102) 

Here U is the universal massless scattering kernel, and 5 is the "massless convolution" in the n variable. 
The solutions are required to obey the asymptotic conditions at n — > oo 

hm flM = V2, l im !?M = i. (103) 

Numerically, the solution of the equations is straightforward in the domain < \ < tt/3 (the sign of \ 1S n °t 
important). In particular, all solutions ei t2 (K) are then real on the positive real axis of k. 

Note that these TBA equations, from a knowledge of the functions on the real line, trivially define e± i2 (k) in 
the domain |Arg(/t)| < 7r/4, which we call the first fundamental domain. This is due to the particular form of the 
integration kernel in the TBA equations, whose poles lie on the aforementioned lines. To obtain the TBA functions 
outside of the first fundamental domain for a given value of \ from a knowledge of their values within the first 
fundamental domain, one needs to perform an analytical continuation procedure (the details of which are to be found 
below). Note that as a function of x, things will get even more complicated. As we will see below, this will require 
another type of analytical continuation, with multiple steps, generalizing the continuation proposed for the case A = 1 
in [9]. 

The first bit of information we can gain comes from looking at the limit k — > of the universal form. Defining 
lim K ^ e i = Li, we get the equations 

Li = iln[l + 2cos3xe L2 + e 2ia ] , L 2 = i In [l + e Ll ] . (104) 

Using the notation Xi — e Li , we can rewrite these as 

X!=xl-l, x 2 {x\ - 3) = 2 cos 3%. (105) 

These are useful equations allowing us to check the numerics. However, these give us the first indication that not 
everything is fine. The third-order equation for x 2 yields three different solutions: x 2 = 2cosx, — cosx ± \/3. For 
Ixl < ""/S, the appropriate one to pick is 

x 2 = 2 cosx- (106) 

At x = ""/S, this therefore means that X\ vanishes, and therefore L\ — > — oo. A naive solution to the TBA equations 
would produce a cusp in x here. If one were to naively solve the TBA for \ beyond this value, one would simply 
regenerate previous results: for x = tt/3 + f], the curves would be identical to those for x = t/3 — r\, which is patently 
incorrect in view of the cusp. The same kind of thing happened for A — 1 at \ — 7r/2,and we refer the reader to the 
discussion in section 3 for a description of the problems and their proper solution. 
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As we expect analytic behaviour as a function of x, some form of analytic continuation will have to be performed. 
Note also that xi vanishes at x = 7r/2: this further problem is not repaired by the surgery at x = 7r/3, so the 
analytical continuation will involve more than one step. 

To make progress on the problem, let us start by reformulating the TBA equations into a more convenient form. 
Introduce the following functions (lu = e J7I 7 4 ) 



P(k) 
Q(k) 



( K -A 3 / 2 )( K -A 3 fi /2 ) 
( K + Af)( K + Af )' 
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{k + wA* /2 )(/s + OjA 3 / 2 ) (k + w Af)(/t + uAf) 



1/2 
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cos6>(k,x) 



A|AfjCos3x 



( K 4 + A|) 1 /2( K 4 + A«)V2- 

The eigenvalues of K\K r are equal to Q(«)e T * e ( K ' x ) . Introducing then the Y functions through 

e ei (n) = P(k)Y!(k), e £2 ( K ) = Q(K)e* i6{K ' x) Y^{ K ) 

the universal form of the TBA can be rewritten 

lnYi(K) = - lnP(/t) + U * (lnQ 2 ( K ) + ln[l + Y+{n')] + ln[l + F 2 " («')]} , 
lnK^/t) = ±i6»(«;,x) - In Q(k) + [/ * {lnP( K ) +ln[l + Yi(k)]}. 

We can in fact verify the following integrals, for d, n e M + : 
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(111) 



where we have to make a branch choice in the last equation. In fact, knowing that the phase of Y 2 is given by ±6, 
we can see that the appropriate choice should yield a simplification of the TBA equations with the identities 



U* lnQ 2 (n) = InP(K) +\uR(k), 
U * ]hP(k) = lnQ(K), 



where 



^_ " 2 + A i« 2 + A « 

H[ >- « 2 -A|« 2 -A 3 fl ' 

The final form of the TBA is therefore 

In Y\ (k) = lnR{n) + U * {ln[l + Y+(k)] + ln[l + F 2 "(k)]} 

2 

The ground state energy then reads 



(112) 



(113) 



lnY^(K) = ±*fl(/s, x) + 17 * In[l + yi(/s)]. 



E ° = 4^rJ d «{^M 1 + 1 i~ 1 (' s )] +Ml + ^2 + '(«)]+ ln[l + K,- '(k)]} 



(114) 



(115) 



Note that, one has to be careful along the lines arg(«) = ±7r/4 where the universal kernel leads to singular integrals, 
which have to be evaluated by contour deformation. First, note that simple manipulations yield 



P(uk)P(cuk) = Q 2 (k), 



Q(uk)Q(ujk) = P(k). 



(116) 



and also R(u>k) = R (uk). Observing that 6{ujk, x) = 0(ojk, x) for real k, and considering the evaluation of the 
TBA functions on these lines, yields the F-system 
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ri(w/s)yi(w/s) = [i + y+(k)][i + y 2 -( k )}, 

Y^(uk)Y^(Cjk) = 1 + Yi(k) 



(117) 



which coincides with the D$ F-system in [29]. In particular, the nontrivial periodicity in the complex plane of k can 
be read from those. Namely, one can show by direct substitution that 

Yt(uj 6 K) = Yi(k), Y±(uj 6 k) = Y±(k). (118) 

This means that the Y- functions naively have Laurent expansions in powers of k 4 / 3 . In particular, the boundary 
pairings appear in this expansion in powers of A| R - Once again, like at the free fermion point, the TBA does not see 
clearly the term proportional to A^A^ in the perturbative expansion of the ground-state energy at weak boundary 
pairing. 

In view of the periodicity of the Y functions, and the consequence this has on the structure of their Laurent 
expansions, we see that the ground-state energy has a naive expansion in powers of A| R . However, as we have stated 
above, problems are encountered at \ = ir/3 which require modifying the TBA in a nontrivial way. We now discuss 
how this surgery is to be performed to ensure analyticity in the boundary phase difference \. 



B. Scenario for the analytic continuation 

As we have argued above, we encounter various types of singularities in the (solutions to the) TBA equations as a 
function of \- Let us define the following terminology: types Ki,K 2 singularities are defined for k obeying 

1 + Yi(k)=0, 1 + Y 2 ± (k) = 0. (119) 

We also define types F± , F 2 ± by 

Y 1 (k)=0, Y±(k) = 0, (120) 

and types G\ , Gf by 

Y 1 -\k) = 0, (Yfr^O. (121) 

Types K\, K 2 and Fi, F 2 have consequences on the ground-state energy integral, whereas G\, G 2 types are harmless. 

We now have to ask ourselves what singularities appear in our equations, and how we have to take them into 
account. From the numerical solution of the TBA equations, we know that at \ = 7r/3, there are singularities of types 
and iq at k = 0. We propose the following scenario for the analytical continuation of the TBA beyond x = i"/3- 

First of all, notice that a singularity of type (as we know appears at \ — 7r 7 '3) requires In — in mod(27r). In 
particular, the real part of lnY 2 must vanish (this is equivalent to saying that e 2 must be purely imaginary, as was 
argued for the case A = 1 at the end of section 3). By examining the form of the universal scattering kernel, and in 
view of the comments made earlier about the use of the TBA equations to write the functions Y in the complex plane 
of k from a knowledge of their value on the real line, we can easily convince ourselves that the condition 3? In Y 2 = 
could only be satisfied on the axes uk, ljk with k£I. We therefore state 

Proposition 1: at \ = 7r/3, singularities of type K 2 and F± occur at k = 0. Increasing \ beyond it/ '3 moves the 
K 2 singularities to wk s , cjk s , with n s e R > 0. The F\ singularity moves to n s (that is, the K 2 singularities are at 
angles ±7r/4 in the complex plane of k, and the F\ singularity sits on the real line). Note that the relative position 
of the singularities is dictated by the F-system. 

If such a singularity structure occurs, the first TBA equation needs to be modified to take the pole contribution 
into account: 

In Y\ (k) = In T(/c, K al ) + Ini2i(/s) + U * {ln[l + Y+(k)] + ln[l + Y 2 {k)]} , (122) 
where we have defined the function 

2 2 

lnT( K , Ks )=ln^ 2 —^| (123) 

as the residue that pops out of the universal kernel. Note that this modification ensures the presence of a F\ 
singularity on the real line at k s , as required. Note also that the extra term does not change the F-system structure, 
since T(ujk, k s )T(luk, k s ) = 1. 

The ground state energy changes under this continuation: we will treat this after solving the singularity problem 
of the TBA equations. 



20 



The major remaining problem is to determine the value of k s . This comes from solving the singularity condition 
obtained by explicitly requiring e.g. Y 2 + (ojk s ) = —1 under Proposition 1. This translates into 



tt (mod 2n) = 9(ujk s ,x) + 



2 - P r 

k Jo 



dn' 



ln[l + Yi («')]. 



(124) 



Note that the fact that Y\ has a type Fi singularity at n s makes the principal part integral well-behaved. The TBA 
system now comprises three equations, i.e. the two basic TBA equations, with the addition of the first singularity 
condition. 

The first thing to notice is that this modification cures our problems at \ = n/3 with the k — > asymptotics. That 
is, the relevant formulas become 



L x = ln(-l) + ^ln (1 + 2cos3 X e i2 + e 2La ) , 



L 2 = -ln(l + e L >) 



Again using the notation x^ 



we can rewrite these as 



Xi 



1 



■ / 2i 



x 2 {x\ — 3) = 2 cos 3x- 



(125) 



(126) 



providing a smooth curve for x± (note the change of sign) around the singular point. 

Now as we tune x U P beyond this, we encounter the possibility of a type K\ singularity, so this is by far not the 
end of the story. From the asymptotics at k — > 0, we expect this to happen at x = 7r /2, where x 2 is bound to change 
sign (meaning that the function Y 2 must become negative on at least part of the real axis of k) . We therefore expect 
a Ki singularity at n — for x = 7r /2- 

Moreover, a careful study of the first singularity condition above shows that one can solve it for k s real only provided 
(here, for A L = A R ) 



Ks < A 3 / 2 



cos ■ 



3x, 



(127) 



The immediate consequence of this is that the window of possible values of the first singular point grows between 
X = 7r/3 and tt/2, but decreases again afterwards. We therefore expect the first singularity to collapse back to zero 
before x reaches n. 

If such a second singularity crosses the real axis, then we have to modify the TBA again. Namely, we now get a 
modified second TBA equation: 



lnY^K) = ±i9(K, x) + lnT( K , n s2 ) + U * ln[l + Y 1 (k)} 
where the parameter k s2 is a solution to the second type of singularity condition 

tt = i lnT(wK s2 , k s ) -ilnR^WK^) - -P / dn' 4 s \ [ln[l + Y+(k')} + ln[l + Y 2 ~ («')]] 

7T J K' - <o 



(128) 



(129) 



The final modification from the second singularity occurs in the first singularity condition: namely, the integral 
kernel appearing in it again diverges, and has to be continued. This produces the same king of additional contribution 
as appeared in the TBA equations, corresponding to the residue of the universal scattering kernel evaluated at 
appropriate points. For clarity, we rewrite the full TBA, including the appropriate branch choices, for the case A = 2. 
In this, we have found it convenient to reabsorb the phases ±i9 and the kernels T in the F-functions. The full TBA 
finally reads 



IilY^k) = In R(k) + U* In 
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A|A3jCos3 X 



(130) 
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(131) 
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with K s i — for x < tt/3 and \ > 27r/3, and n S 2 — for x < 7r /2- The evolution of the singularities is as follows: 
k s i first appears at x = 7r/3, and moves up. At x — 7r /2, k S 2 makes its appearance, and starts moving up. As it 
does so, it pulls k, s \ back down, until the latter vanishes again at x — 27r/3. After that point, only k S 2 remains. This 
behaviour is illustrated in figures (4, 5), where the two singularities are plotted as a function of the phase difference 
X for various values of the boundary parameters A L = Ar. 




FIG. 4. The singularity k s i plotted as a function of the phase difference \ f° r seven values of the boundary parameters 
A L = A R . 




FIG. 5. The singularity k 3 i plotted as a function of the phase difference \ f° r seven values of the boundary parameters 
Ah = A R . 



The ground state energy is given in these notations by the integral expression (115) analytically continued for the 
singularities of type and Fi at x = an d those of type K\ and F 2 ± at x = 7r /2- The extra kernels from 
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these respective divergences cancel pairwise, and one is left after a few basic manipulations with the modified integral 
expression for the ground state energy (valid for all values of \) 



En — 



-I 

4^R 



+ V21n 



dn { In 



1 + 
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1 Y 2 (k) + 



v s2 
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\sl 



- V21nYi( K )|. 



(132) 



Due to the fact that we have reabsorbed the singularity contributions into the Y functions, they are always real and 
positive on the real line. 

The numerical solution to the full TBA equations (including the singularity conditions) over the full interval of \ 
from to 7r turns out to be a considerable challenge. The evaluation of the principal part integrals requires many 
sampling points, but this isn't the main problem. The main problem is that the numerics tend to be unstable within 
a small domain around x ~ 1-9- This is due to the fact that the right-hand side of the singularity conditions arc 
nonmonotonous functions, from which it is therefore difficult to isolate the correct root. Four integral equations have 
to be simultaneously obeyed, so the iteration procedure in general has many unstable directions. A simple root finding 
algorithm is insufficient, and one has to make the code rather elaborate to obtain sensible curves. 

The curves for the ground-state energy as a function of \ f° r various values of the boundary parameters A are 
presented in figure (6). ^From these, it is transparent that our construction provides a completely interpolating solution 
between the ~ cosx behaviour at small A to the <~ x 2 one for large A, as expected from the appropriate conformal 
limits. For large A, the fit is perfect with the expected functional form (easily derived from mode expansions) 



7T 3x£ 
24i? + 4tt 



(133) 



with the last equality valid for the present case A = 2. 
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FIG. 6. Ground-state energy of the two-boundary sine-Gordon model plotted as a function of the phase difference x f° r 
seven values of the boundary parameters Al = Ar. The interpolation between the conformal limits is clearly seen. 



C. Some remarks on the general singularity structure 



We are not sure what the general singularity structure will be for other values of A. A generalization of the Ising 
model argument given in section II indicates that there always will be a singularity at ^. Indeed, following [12], 
boundary states are properly understood by considering a semi-infinite system, and choosing at infinity cither the 
boundary condition <f> = or <p = ^ (the equivalent of spin up or down in the Ising model). For \ e [0, §] , the 
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first condition defines the true ground state (equivalent to h > in the Ising model) while for x e [ , 7r] , the first 
condition corresponds to an excited state. Equivalently, if one sets always <j) — at infinity, then the case x£ [f , 7r] 
should correspond to an excited boundary state. 

We pause here to stress that the ground state energy does not have any singularity, and that there is no level 
crossing the ground state at any of the singularity values in the massless, two boundary problem. Indeed, the 
foregoing argument assumed values at infinity fixed by the bulk interaction term, with periodicity <f> — > <f> + jj.. In 

the massless two boundary problem, the periodicity is <f> — > <j> + so the first level crossing is expected to occur at 
X = 

The existence of the singularities we are struggling with translates into interesting behaviours for the excited levels 
of the hamiltonian. This can be seen quite clearly in the case of the Ising model, for instance with boundary fields 
hi = h fixed, h r = h' varying. Since the ground state is analytic, the first gap behaves as at small magnetic 

fields, and pinches the ground state in a V shape around the symmetric point h' = 0. 

Technically, the existence of singularities in the TBA approach follows from the dependence of the scattering matrices 
on (A + l)x instead of \- This dependence is a consequence of the fact that the scattering matrices are in essence IR 
defined. Starting from the UV action (1) the Dirichlet IR boundary conditions are approached along the operators 

"dual" to e 1 ^ 2 , i.e. e lfid< t> d / 2 where cj> d = cf> - 4> and |& = || . It follows that the dual coupling A d cx (a) " and 

therefore, if x is the phase of the UV coupling, the phase of the IR coupling becomes (A + l)x- Clearly, the IR action 
(or the scattering matrix) by itself is not sufficient to describe all values of x Mod 2tt, and it is natural to expect 
that the reduced information contained in (A + l)x Mod 2-7T has to be supplemented with specification of a boundary 
state. While it is possible to associate a conformal Dirichlet boundary state with any value of the field (f)(0), it is well 
known (see eg [27], whose conventions we follow) that for a given radius R = j|, there are, when R 2 = 4j£r, A + 1 
"special" boundary states. Presumably, the IR action (or the scattering matrix) supplemented by the choice of one 
of these boundary states would allow to explore the whole domain \ Mod 27r, which would correspond in the TBA 
approach to analytic continuation at the values x = jzp[- 

It may also be that more values require continuation. This can be inferred from the limit k — > 0, which gives the 
following values of the Xj for general A: Xj = ■ We thus expect from this to observe singularities also at 

values -^-j- , f , § , and multiples (the pattern is similar but different from the one of excited boundary states discussed 
in [6]). More work is required to clarify this question. 

VI. THE DOUBLE KONDO MODEL 

The TBA approach to the double Kondo model is bound to be rather intricate, as the boundary scattering involves 
now (except for spin 1/2) boundary degrees of freedom, whose inclusion into the boundary state and the subsequent 
R-channel TBA is not obvious. Fortunately, the problem is perfectly well suited, on the other hand, to the Destri de 
Vega approach. 

The general framework to apply the DDV approach to a problem with two boundaries has been discussed in 
details in [2]. There, the starting point was the inhomogeneous 6 vertex model with boundary fields. The role of the 
inhomogeneitics was to introduce a bulk mass scale, while the boundary fields constrained the value of the sine-Gordon 
field on the boundaries. The net result was an expression for the ground state energy of the sine-Gordon model with 
double Dirichlet boundary conditions. 

It is well known how to introduce a different kind of inhomogeneities in the 6 vertex model to get a theory that is 
massless, with free L movers and R movers interacting with a Kondo type impurity. It is similarly possible to obtain a 
theory with free L movers and R movers interacting with two different Kondo type impurities. Forgetting about the L 
movers, one can then transform this problem by folding into the problem we are interested in. It is a straightforward 
calculation to then obtain expressions for the ground state energy - these were actually almost already written in 
[2] 2 . Since the most interesting limit is the isotropic one (XXX model), we give results in the repulsive regime of the 
sine-Gordon model, parametrizing fi 2 = 8tt^-. One finds then 



1 - e 



f(0'+iO) 



/OO 
d0'$(0 - 6>')Im In 
-OO 

9 f°° r i 

E=-—\ p^Imln l-e^^ ' (134) 

J-oo L J 



2 There is a misprint in this paper, as equation (9.3) corresponds to Neumann, not Dirichlet boundary conditions on the other 
boundary. 



24 



where $ is the kernel already encountered in section IV 



m = r *L f nHt - 2) \, e 2 ^ (135) 

W J_ 00 2tt 2 coshxsinh(t- l)x V ; 

and the boundary kernel is 

= l L V 2cosh,sinhft-i), sin2 ^ - - y 7-r + ^ r ) ( 136 ) 

The phase is such that at small coupling 9 B — > — oo the boundary term vanishes while at large coupling 9 — > oo it 
reads iPtdry = n + T^)' this ma ybe modulo 27r. 

Note the deep similarity of these equations with the ones usually written in the bulk - in the latter case, there 
would be a ioj term instead of a iPbdry term, the u> corresponding to a soliton fugacity. 

While bulk equations involve a constant ui, the variation of Ptdry with the rapidity leads to some obvious difficulties. 
Let us consider for instance the simplest case of a spin 1/2 Kondo impurity at one boundary, Neumann conditions 
at the other boundary, and restrict also to the isotropic case t = oo. As discussed in [8] going from UV to IR in this 
case amounts to fusing the identity (Kac Moody) representation that initially goes through the system with the spin 
1/2 (Kac Moody) representation, so the ground state energy in the IR should correspond to the conformal weight 
h = l/2(l/2+l)/l + 2= 1/4. This is in agreement with our formula, as the usual dilogarithm analysis gives 

^-\{^)'-\ d") 

If however we put two Kondo impurities, we see that our formula gives hm = 1, while fusing with the Kac Moody 
twice gives the Kac Moody identity representation back, with lowest conformal weight h = 0. The h = 1 weight is 
instead an excited weight, and therefore our DDV equation, while it captures the ground state at small coupling, 
captures instead an excited state at large coupling. 

It was also argued by Affleck and Ludwig in [26] that the picture for j = \ essentially carries over to the case of 
arbitrary j: in cither case, flow from UV to IR is described by the absorption of an eletron by the impurity, and a 
phase shift of ir/2 for the fermion wave function. Hence all the foregoing remarks extend to the case of higher values 
of the spin as well, which agrees with our general formulas, as the value of Pbdry(oo) does not depend on the spins 
in the isotropic limit t — > oo. In terms of the free boson, the boundary condition corresponding to the Kondo strong 
coupling fixed point is <pi — <fi r = -^/l^niodulo \/2tt. 

What happens in the anisotropic case is somewhat less clear, as the boundary condition left over in the strong 
coupling limit now seems to depend on the impurity sl(2) q spin. 



VII. CONCLUSION 



It is tempting to carry out a little bit further and propose a DDV equation for the double boundary sine-Gordon 
model. 

To do so, we would like first to get back to the free fermion case. We introduce the two eigenvalues of the matrix 
KiK ri which we call A±. Their expressions are (setting fi = 1) 



Ah 



1 



e 2f> + AfA 2 r cos 2 X ± i\J AfA 2 sin 2 2 X ~ e 2e (Af + A 2 r + 2AfA 2 r cos 2 X ) 



{e e + Af)(e $ + A 2 .) 
and we rewrite the ground state energy (60) as (recall k = e e ) 

E = -—J d9e e \n(l + A + {9)e- Re )+-—J d9e e In (l + A-(9) e - Re ) 



(138) 



(139) 



We now perform a shift of the variable of integration, 9 — > 9 + ^ — it in the first integral, 9 — > 9 — ^ + it in the 
second integral. Here, the t's are necessary to avoid the zeroes of the logarithms in the integral. Care must also be 
exercised in moving the contours because of the branch cuts in the definitions of the square roots in A± . This gives 
rise to the new result 



E = 



1 

4l7T 



e% 1 



1 

AT 



J C*2 



(140) 
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where we have introduced the eigenvalues of the product RiR r (no complex conjugation, see later), 
1 



(e« - iAf)(e e - tA 2 ) 
and the contours are represented in the figure 7. 



e 26 + A 2 A 2 cos 2 X ± ^ V A " A " sin2 2 * + e2 "( A l + A r + 2A ' A r cos 2 X ) 



(141) 



FIG. 7. Contours used in equation (141). 

Similar manipulations could lead to the same result with A + replaced with A_ , and the two are in fact equal. After 
some simple redefinitions we thus obtain the result 



E : 



4j7T 



f e*ln(l-e'W) + f e e ln(l-e-W) 



where 



f{6) =Rie 8 -lnA± 



(142) 



(143) 



We conjecture that the DDV expression in the general case is very similar, with the exact same form for the energy, 
and as for the source equation 



/(*) 



Rie e - In A+ + f $(0 - 9') In (l - eJ) + f &(6-#) In (l - e" / ) 



(144) 



and <E> is the kernel already used in the DDV equations of the previous paragraph. 

Here, A± are eigenvalues of RiR ri which have a structure entirely similar to the one for the case A = 1, 



A± 



1 



e 2X8 _ A A+1 A A+1 CQS71 ± . [ A 2A+2 A 2A+2 ^2 J] + £ 2X6 ^X+2 ^X+2 + 2 COS 7? A A+1 A A+1 )] V * 

(and recall r\ = (A + l)x)- An important property is that \+(0 + in) = A_(#). 
It is interesting to see what our conjecture gives in the case where 1 — > 87r. There one finds 



(145) 



A± 



-1 - A ( A r cosx± i\J AfA 2 . sin 2 x + A 2 + A 2 - 2A ( A r cos rj 
v/(l + A 2 )(l + A2) 



(146) 



In this limit, the dependence on the energy of the extra soliton fugacity has disappeared, and the ground state energy 
follows from the usual dilogarithm calculation 



E 



24R 



i 6 ' 

l — or 



The eigenvalues being pure phases, one has 



a = ±arg A± = =p arctan 



A 2 A 2 sin 2 x + A 2 + A 2 - 2A; A r cos x 



1 + A;A r cosx 



(147) 



(148) 
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The relation between the coupling A and the bare coupling is simply A = ttA. 

We check this expression by comparing the resulting "current" in applications to Josephson juncionts to the one 
obtained in [25]. The current follows 

da . . a ,...„.. 

Joe a— oc A; A r sin x (149) 



dx 'i&i&l sin 2 X + A 2 + A 2 . - 2A ; A r cos x 



and a can be written 



1 + A;A r cosx 



a = arccos ^= | = == (150) 



It is easy to check that this coincides with the formula in [25] but for the redefinition 

a ^there 

&here = _ T 2 (151) 

1 ^there 

Getting back to arbitrary values of A , in the limit J? — ► oo, it is easy to check from the same kind of calculations 
that the correct value 

^—M 1 - 6 ^) (152) 

is recovered. 

Of course, these few tests are far from enough to completely justify our proposal, and we expect to get back to it 
in more details in the future. 

As we were completing this paper, a preprint appeared [31] where the same problem is discussed, and the problem of 
analytic continuation is also claimed to be solved. We believe that the solution proposed in [31] is actually incorrect, 
for technical reasons we now discuss. 

In this preprint, the authors claim to find the solutions to the singularity conditions on the imaginary axis of k, for 
any value of A. While this is certainly true for A = 1 as shown in our earlier work [9], we have shown that for A = 2, 
the singularities should be found on the lines with angles ±7r/4 in the complex plane of k, instead of the imaginary 
axis. Moreover, the authors of [31] argue that only one singularity shows up for any value of A, whereas we have shown 
(using for example of k — > asymptotics and the V-system arguments) that two singularities have to be treated for 
A = 2 to ensure analyticity in x- As an illustration, consider equation (3-12) of [31], which is the universal form of 
the TBA corrected for the first singularity. The result for the correction is the residue InS'ao, whereas we correct 
with the residue of the universal scattering kernel (see our equation (122)). The reason for the difference seems to be 
that the authors of [31] have not taken into account the fact that the function t\ diverges at k = when x = t/3, 
in addition to the first logarithmic kernel involving e 2 (the notations are different: eo in their work is what we have 
called £2 here). These two divergences each modify the TBA equations, and correcting for only one of them leads to 
the wrong results. This is transparent if the whole analytical continuation procedure is done using the universal for 
the the TBA equations from the start, as we have done. 



ACKNOWLEDGMENTS 



We thank P. Dorey, R. Egger and A. Zamolodchikov for helpful conversations. This work was supported by Stichting 
FOM, the DOE and the Humboldt Foundation (HS). 



APPENDIX A: THE FREE FERMION POINT 



At /3 2 = 4?r, the sine-Gordon action has as is well-known a direct correspondence to the action of free fermions. 
In order to correctly treat our double-boundary theory at this free fermion point, we have to be very careful with 
our bosonization/refermionization identities. These are much simpler in the bulk, where we can simply toss away the 
finite-size terms, and for simple periodic boundary conditions, for which the quantizations rules of the zero modes can 
be directly read out. We will start from the very basics to illustrate the procedure forced upon us by the two-boundary 
geometry we are considering, and give every detail of the derivation from the canonical quantization of the free boson 
to the final formula for the ground-state energy. 

Consider thus the interval J:i£ [0, R] on the real line. We want to study a massless bosonic field <p(x) which is 
free in the bulk (i.e. for x G ]0, R[) but has some boundary contribution to the action at the points x = 0, R. If we 
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intend to develop <fi(x) in a mode expansion, we might very well be tempted to use periodicity with period R. As we 
do not wish to identify 4>{x) at the points and R, in order to accommodate different boundary effects on both sides, 
this is not flexible enough. Instead, we will extend the definition of cf){x) to the interval x € [0, 2R] by using 

tj>{2R-x) = <f>{x), xe[0 7 R}. (Al) 

We can then consistently extend the definition of <p(x) to the entire real line by using 

4>{x + 2R) = 4>{x), xeR. (A2) 

This is important: the requirement that the fields on both ends of the original interval were not identified with one 
another, has required us to define the field as periodic on an interval with double the original length. 
Let us consider the free Hamiltonian 

I r R 

H °=2j +(M) 2 ]- (A3) 

Imposing the properties of 4> in (A1,A2), we can write the mode expansion 

4>(x, t) = 0° + n A + * J2 -an cos ^ e -^l R (A4) 

where the commutation relations of the modes read 

[0 o ,ITo] = i, [a n ,a m ] = n8 n +m,a (A5) 

with all others vanishing. In terms of these, the Hamiltonian reads 

tfo = ^| + ^£ a « a — (A6) 

Alternately, we can play another game and define a chiral left-moving boson as 

Mt + x) = <t? L + n ^±f + -i= E ~a n e—^/ R (A7) 

obeying the quasi-periodicity relation 

4> L (x + 2R) = cl) L {x)+n . (A8) 

We can without contradiction do a nonlocal identification between our original bosonic field and this newly invented 
left-mover. Namely, we may at leisure impose the operator identity 

<t>(x, t) =4> L (t + x) + <t> L (t - x) + 4>° - 2<jy° L (A9) 

where the periodicity requirements of the original boson are automatically fulfilled given the quasi-periodicity of the 
chiral one. Note that under these circumstances, the Hamiltonian (A6) becomes 

Ho = / dx{d x <t> L ) 2 . (A10) 
Jo 

Moreover, it is a trivial exercise to show that taking 0° = 4^/2 leads to consistent canonical commutation relations 
for the chiral field: 

[H{x), d x ,<f> L (x')} = l -J2S(x-x' + 2Rn). (All) 

Let us now move on to the model we are interested in. At the free fermion point, the boundary contributions to 
the Hamiltonian read 

H B = Ai cos V^(f>{0, t) + A r cos(V^>(-R, *)-*)• ( A12 ) 

To refermionize this at 1 — Air, we start by noting that the original boson at x = 0, R is simply expressed in terms 
of the chiral boson: 
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0(0, t) =2<j> L (t), 



<j>(R,t)=2^ L (R + t)-U Q 



(A13) 



Progress is now simply a matter of using the definition of the chiral fermion from our favourite bosoniza- 
tion/rcfcrmionization rules, 



* = 



f 



^Ait a 



/4w<t>L 



to obtain the equivalent fermionic boundary Hamiltonian 



H B = ^=[*(°) + * f (°)] + ^[f*(RY x + fW(R)e- ix ] 



(A14) 



(A15) 



where we have defined / = e ~ l ^ Uo and rescaled our boundary parameters as A = y/n A, setting the cut-off a equal 
to one. Note that we now face a linear fermionic term at x = 0, and that the behaviour of the operator IIo, and thus 
that of /, is still unspecified. To resolve this is simply a matter of noting that the operator n becomes quantized 
through the periodicity requirement that we imposed on the original theory. Writing the boundary contributions in 
terms of the left-mover, we get the IIo quantization relation and (being very careful with noncommuting operators) 
its fermionic implication 



n = x/ttti, V(x + 2R) = -V{x) 

with neZ. These then tell us that the operator / in fact obeys the Majorana relations 

f = P, {/,/} = 2, {/,*} = 0. 



(A16) 



(A17) 



For aesthetic reasons, we would like our Hamiltonian to be quadratic in fcrmions, and somehow more symmetric with 
respect to the left and right components. In order to do this, we perform the simple transformation 



/ — > i a b 



(A18) 



where a, b are two new Majorana fermions, living respectively on the left and right boundaries. One can check that 
all anticommutators are preserved under this transformation. 
We thus finally obtain a better-shaped real-time action, 



dt 



2R 



dx 



*U(dt - d x )* + - d x )& 



dt 



l -{ad t a + bd t b) + ^o[*(0) - * f (0)] + ^=fe[*(i?)e lx - ^(R)e~ ix } 



(A19) 



Our theory has thus boiled down to something rather obviously tractable, namely one that is quadratic in fermions, 
and which we can consequently solve exactly. Varying the action and eliminating the boundary fermions, we obtain 
two sets of boundary conditions at x — and R: 



*(0+,t) + tf^O+t) = tf((T,t) + * f ((T,t) 



(A20) 



$*(o+t) -d t ¥(o+,t)- d t *(o-,t) + d t &(o-,t) = 
= 4r - * f (o + ,t) + *(o-,t) - **(<)-,*)] 

$(R + ,t) +e- 2tx ¥(R+,t) = «f(R-,t)+e- 2lx ^(R-,t) 



(A21) 
(A22) 



d t V(R+,t) - e- 2ix d t &(R+,t) - dt*(R-,t) + e- 2ix d t &(R-,t) = 
A 2 

= ~f [*(# + , t) - e- 2tx ¥{R+,t) + *(iT , t) - e- 2ix &(R-,t) 



(A23) 



The fermion is a free left-mover in the bulk, with possible discontinuities at x = 0, R. We choose the general mode 
expansions 



k 



< x < R 
R < x < 2R 



(A24) 
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and substitute them back into (A23), keeping in mind the antiperiodicity of the fermions under x — > x + 2R. This 
scheme is consistent provided the momenta k are quantized according to 

1 + W = 0, 

o [4(ifc) 2 + A^cos2 X ] [2ik-Af][2ik-Al] iikR 

[2ik + Af][2ik + A*} + [2ik + Af][2ik + A*} { ' 

The ground-state energy is then given by 

i i r dw_ -, ,■ 

E o = t:Y k = -r- \ dkk-^— = -— / dk\n [1 + W]. (A26) 

where the anti-clockwise contour C surrounds all the roots of our quantization condition (A25) lying on the negative 
half-line k < 0. This is illustrated in figure 8. 



Im k 



X X X X X X X X X i 



■> 



X X X X X X X X 



Re k 



FIG. 8. Contour of integration for performing the summation over occupied energy states when calculating the ground-state 
energy. The crosses represent allowed values of the quasimomentum, which are symmetrically distributed around k = 0. 

The contour C can be split in two different parts, Ci t 2, corresponding respectively to the upper and lower branches. 
The first integral (over C\) then converges after performing a Wick rotation k = ik, whereas the second integral 



(over C 2 ) converges after extracting a multiplicative factor i^rr^jjl 



ik-K\ c -UkR 



pfe+Af]pfc+Aj] 



from within the argument of the 



logarithm while performing the Wick rotation k 
the final formula 



En — — - 



1 

47T 



dn In 



-ik. Setting the nonuniversal term J Q dkk to zero, we then obtain 
2 + A2A2 C os2 x _ Kfl , (K-Af)( K -Aj) 



-2kR 



(K + Af)(K + At) 



(K + Af)( K + Ai) 



+ ei + e r 



where e; jr are the boundary intensive energies 



■— r 

4Wo 



dk 



(A27) 



(A28) 



We have thus obtained the ground-state energy of the (3 2 = 4ir sine-Gordon model in a finite width R, in the presence 
of two integrable boundary interactions. Although the integral is not expressible in terms of elementary functions, it 
gives us the dependence of E in terms of the boundary parameters A;, A r and phase \- 



APPENDIX B: BULK SINE-GORDON SCATTERING MATRICES 



In the bulk, the scattering of the soliton and antisoliton is described by the general processes 

= a(0i - 2 )4(02)4(0i), 

4(0i)4(0 2 ) = a((?i - 02)4(02)4(9,), 

4(0^4(02) = b(0, - 2 )4(0 2 )4^i) + <0i - 02)4(^)4(h)- (bi) 

The amplitudes for these scattering processes are well-known [21]. When 1 = -£fj, A e N, the theory contains the 
following fundamental particles: the soliton, the antisoliton, and A — 1 breathers of mass m n — 2msm^. The bulk 
scattering simplifies since the backscattering amplitude c(9) vanishes; the other two amplitudes read 
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a{9) = isinh[\{8 - m)]p{8), 
b(8) = -ismh[X9]p(9), 

<fi\ ^mm-u^TYi x .X9 f* F l {8)F l {m-9) 



r(2ZA + i^)r(l + 2ZA + i^) 



r((2/ + 1)A + if)T{l + (21 - 1)A + iM) 
A useful integral representation for this is [23] 

dy 2A0y sinh[(l - A)y] 
sinhy cosh Ay 



<a\ 1 (■ r d y ■ 2X6 

^ = sin[A( 7 r + ^)] eX HV- 00 % Sm - 



(B2) 



(B3) 



At the reflectionless points 1 = , A G N, the bulk scattering amplitudes simplify considerably. Namely, the 
backscattering amplitude c(8) vanishes identically, while the others simplify according to 



^)=-ri- os(s+z 



j=i cos - 4) 



(B4) 



as 



sin(A(7r + i6')) 

These have to be supplemented with the scattering amplitudes for processes involving breathers, which are defined 
BU9i)Bl{0 2 ) = S^ m \8 x 9 2 )Bl{8 2 )B\ l {8 1 ), 



S {n) (9) 
S {n ' m \8) 



sinh 9 + i cos ff sin 2 (gjgg 
sinhfl-jcosff 11 sin 2 ( (!i_^_|_,|)^ 

sinhfl + i sin(^±^) sinhtf + i sin(^^) 
sinhfl - zsm(^^) sinlifl- isin(i2^L) 

-n— 
~2\ 



X 



2A y " " oiii V 2A 



- sin 2( ( m -»-2/K _ z | )cos2( 



2 f (m+n— 2()7T 
' 1 _ 2A 



-if) 



n> m. 



(B5) 



The integral representations of these amplitudes read 
S (n) (8) =cxp 



ciy . ,2X8y. sinh ny cosh y 

siiif ) 

2y tt cosh Ay sinhy 



S {n ' m) {6) =cxp (i 1°° - 

V ./-oo 2/ 



dy . 2X9y 



sin( 



&n.m 2 



cosh y cosh[(A — n)y] sinh my 



cosh Ay sinhy 



n > m = 1, A — 1. 



(B6) 



APPENDIX C: BOUNDARY SINE-GORDON SCATTERING MATRICES 



1. Solution to the boundary bootstrap 

The amplitudes for the scattering of the soliton and antisoliton on the boundaries were computed in [12]. The 
bootstrap solution yields these amplitudes in terms of two parameters 77, 7? whose relation to the parameters in the 
original Lagrangian was not specified in [12]. For now, we state the results in terms of 77,1?. For general A, the 
amplitudes for soliton and antisoliton boundary scattering read 

P±{8) — (cos(— i\9) cos 77 cosh $ =p sin(— i\9) sin 77 sinh 7?) R(8), 

Q±{9) = ~ sm(-2iX9)R{9), (CI) 

where 

R(9) = R (9)R 1 (9). (C2) 
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The first of these functions is given by 

F o (0) 



Ro(0) 



F (-ey 

where 

r(4fcA + i2M) r(l + 4fcA + z^)r((4fc+l)A)r(l + (4A;-l)A) 
Fo( ) ~ r(A + i2M) AA r((4fc + 1)A + T(l + (4fc - 1)A + i^)r(l + AkX)T(4kX) ' 

Its integral representation reads 

dy . 2X9y sinh((A - l)y/2) sinh(3Ay/2) \ 



J?oW=exp f-i r^sin(^) 

V i-co 2/ 



The second function reads 
where 

a(x,6) = 



Ri(9) 



sinh(y/2) sinh(2Ay) / ' 

a(r),9) a(i$,6) 
cos r\ cosh d 



u(x, -e + if )ii(-x, -0 + if )n(x, 6» - if )n(-s, - if) 



n 2 (x,if)n 2 (-a; 

n , T m fr £(| + (gj + |)A + I + ig)r(j + (2/ + |)A + |) 
1 ' j f = 1 r(i + (2Z + |)A + | + zM)r(I + (2; + I)A+f)' 



This obeys the property 

a(x, 9)a(x, —9) 



cos(x + i\6) cos(x — iX9) 
(correcting the equation in [12]). The integral representation for these is 

, m ( f°° ^ , i^u ,A0u N cosh(2£M) 

ctO, 61) = exp i / — sinh(Ay(l + — ))sin( — )- 



y it tt sinh y cosh(Ay) / 

Again, at the reflectionless points 1 — ^j,AN, these simplify considerably: 



A-l 



Row = n 



cos(f + 



,=1 + fj[-<f)' 

^ m - cosx A cos - ix - 1 + »f ) 
(,) cos (-^)M cos (^-t-l-^)' 

In the massless limit, we get 

i?0(^-e)=e-^(A-l) 

~ g ) = 2e -A0-*„A/2 

cos ?7 

<7(il?,if -0) 2 -A+l e -Afl/2- 1 J/2 



cosh TT A coshf- - — - i— + i^. - 4Z) 

llj=l^ ot > n \2 2A 4 4A +t 2A 4 4 ■> 
4 e -i^(A-l)-A0-if A 



i?(if -0) 



]Tv =1 e e + ex+^-'TT+'f 
The breathers also scatter off the boundary, according to 

BUe)B = R%\e)Bt(-6)B. 
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Since the breathers are bound states of (multiple) solitons and antisolitons, their boundary scattering amplitudes can 
be worked out systematically. This procedure was carried out in [24], with the result 

R^\e) = R { ; i \e)R[ n) (9), (ci3) 

where 

R (n )(e) = cos(|f - if) cos(|f + f+ i|) sin(f - if) sin(^ - ,0) cos 2 (|| j-f+if) 
cos(^+i|)co8(^ + f-if)sin(f+if)fJsin(^+i(9)cos2(^ + f-if) 

and the function R] n '(9) depends on the parity of n, i.e. 

R (2n) {9) = S ( - 2n \r 1 ,6)S ( - 2n \i#,6), n= 1,2,... < A/2, 

™ cosff - itilZL) + i s inh(9cos(f + itilZL) + i s inh<9 

S (2n) (x,6)=T[ — ^ T . \ . ; ^ Fr 4r J (C15) 

r=i cos(f - ii^) -isinh6icos(f + ^^) - i sinh6» 



or 



Rf n ^\e) = S (2n - 1) (r 1 ,9)S {2n ~ 1 '>(itf,9), n=l,2,... < (A + l)/2, 

S (2n-i)( 6») _ cos(f) +^sinh6> ^ cos(f - g) +zsinh g cos(f + g ) + zsinh6> 
lX ' j cos(f) -isinh<9 11 cos(f - f ) - isinhfl cos(f + f )-isinh6>' 

The unified integral representation for this is worked out to be 



(C16) 



S in >(x,e)=cxp\iTr-i — sin( -) — ; • (C17) 

y 7r cosh \y sinh j/ / 
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2. Parameter correspondence 

A crucial ingredient missing in the formulas above is the exact relationship between the boundary parameters £, k 
coming out of the bootstrap, and the physical variables Ai yr ,x appearing in the boundary Hamiltonian we started 
from. This correspondence was worked out by A. Zamolodchikov [28] (see also the semiclassical treatment in [30]). 
For a problem where the bulk sine-Gordon action has a term 2/xcos/30 and a boundary term 2\ib cos/30/2, one has 



02 

cosh( — (■& ± ir])) = M sin 

8tt yfji 



0? 



(C18) 



In the massless limit, we will in general thus get 



0« (A + l)ln 



2 , / sin ■ 



A + l 



77 = (A + 



(C19) 



Now from the paper [18] we have that, in the massles limit, = 9' — 6b so Tb = \e 0B = me 6 ^ (in the last formula, 
the factor \ comes from your normalizations in the massless limit, different from the ones in [18]). Here m is the 
soliton mass, related with the bare sine-Gordon parameter \i by [16] 



r(i/A + i) 

7rr(A/A + l) 



m^TTT{\/2 + 1/2A) 



2r(l/2A) 



2A/A+1 



(C20) 



Putting everything together just reproduces the relation between Tb and \ib mentioned in the text. 



[1] J. Cardy, "Conformal invariance and surface critical behavior", in Conformed Invariance and Applications to Statistical 
Mechanics, eds. C. Itzykson, H. Saleur and J. B. Zuber (World Scientific, 1988). 



33 



[2] A. LeClair, G. Mussardo, H. Saleur and S. Skorik, Nucl. Phys. B 453 (1995), 581. 

[3] H. Itoyama, T. Oota, "Normalization of Off-shell boundary state, g-function and Zeta function regularization" , hep- 
th/0206123. 

[4] P. Dorey et al, Nucl. Phys. B 525, 641 (1998). 

[5] G. Feverati, P. A. Pearce and F. Ravanini, Phys. Lett. B534 (2002) 216. 
[6] Z. Bajnok, L. Palla and G. Takacs, Nucl. Phys. B622 (2002) 565. 
[7] P. Mosconi, G. Mussardo and V. Riva, Nucl. Phys. B621 (2002) 571. 
[8] I. Affleck, Nuc. Phys. B336 (1990) 517. 

[9] J.-S. Caux, H. Saleur and F. Siano, Phys. Rev. Lett. 88 (2002) 106402. 
[10] K. Bardakci and A. Konechny, "Tachyon instability and Kondo type models", hep-th/0009214. 
[11] J. A. Harvey, D. Kutasov and E. J. Martinec, "On the relevance of tachyons", hep-th/0003101. 
[12] S. Ghoshal and A. B. Zamolodchikov, Int. Jour. Mod. Phys. A 9 (1994), 3841. 

[13] P. Baseilhac and K. Koizumi, "Sine-Gordon quantum field theory on the half-line with quantum boundary degrees of 

freedom", hep-th/0208005. 
[14] C. Destri and H. de Vega, Phys. Rev. Lett. 69 (1992) 2313. 

[15] P. A. Pearce and A. Kliimper, Phys. Rev. Lett. 68 (1991) 974; A. Khimper and P. A. Pearce, J. Stat. Phys. 64 (1991) 13. 
[16] Al. Zamolodchikov, J. Mod. Phys. A10 (1995) 1125. 
[17] R. Chatterjee, Nucl.Phys. B468 (1996) 439. 

[18] P. Fendley, H. Saleur and N. Warner, Nucl. Phys. B430 (1994) 577. 
[19] H. Au. Yand and M. Fisher, Phys. Rev. Bll, 3469 (1975). 

[20] L. Takhtadjian and L. Faddcev, Theor. Math. Phys. 21, 160 (1974); V. Korepin and L. Faddeev, Theor. Math. Phys. 25, 
147 (1975). 

[21] A. B. Zamolodchikov and Al. B. Zamolodchikov, Ann. Phys. 120, 253 (1979). 
[22] M. Ameduri, R. Konik and A. LeClair, Phys. Lett. B 354 (1995), 376. 
[23] P. Fendley, H. Saleur and N. P. Warner, Nucl. Phys. B 430, 577 (1994). 
[24] S. Ghoshal, Int. Jour. Mod. Phys. A 9, 4801 (1994). 

[25] I. Affleck, J.-S. Caux and A. M. Zagoskin, Phys. Rev. B 62 (2000), 1433. 

[26] I. Affleck and A. W. W. Ludwig, Nucl. Phys. B352 (1991) 849. 

[27] I. Affleck, M. Oshikawa and H. Saleur, Nucl.Phys. B594 (2001) 535-606. 

[28] Al. B. Zamolodchikov, private communication. 

[29] Al. B. Zamolodchikov, Physics Letters B 253, 391 (1991). 

[30] A. Chenaghlou and E. Corrigan, Int. Jour. Mod. Phys. A 15, 4417 (2000). 

[31] T. Lee and C. Rim, "Thermodynamic Bethe ansatz for boundary sine-Gordon model", hep-th/0301075. 



34 



